THE RANDOM-CLUSTER MODEL 



Geoffrey Grimmett 

Abstract. The class of random-cluster models is a unification of a variety of sto- 
chastic processes of significance for probability and statistical physics, including per- 
colation, Ising, and Potts models; in addition, their study has impact on the theory of 
certain random combinatorial structures, and of electrical networks. Much (but not 
all) of the physical theory of Ising/Potts models is best implemented in the context 
of the random-cluster representation. This systematic summary of random-cluster 
models includes accounts of the fundamental methods and inequalities, the unique- 
ness and specification of infinite-volume measures, the existence and nature of the 
phase transition, and the structure of the subcritical and supercritical phases. The 
theory for two-dimensional lattices is better developed than for three and more di- 
mensions. There is a rich collection of open problems, including some of substantial 
significance for the general area of disordered systems, and these are highlighted when 
encountered. Amongst the major open questions, there is the problem of ascertaining 
the exact nature of the phase transition for general values of the cluster-weighting 
factor q, and the problem of proving that the critical random-cluster model in two 
dimensions, with 1 < <? < 4, converges when re-scaled to a stochastic Lowner evolu- 
tion (SLE). Overall the emphasis is upon the random-cluster model for its own sake, 
rather than upon its applications to Ising and Potts systems. 
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1. Introduction 

During a classical period, probabilists studied the behaviour of independent random 
variables. The emergent theory is rich, and is linked through theory and application 
to areas of pure/applied mathematics and to other sciences. It is however unable 
to answer important questions from a variety of sources concerning large families 
of dependent random variables. Dependence comes in many forms, and one of 
the targets of modern probability theory has been to derive robust techniques for 
studying it. The voice of statistical physics has been especially loud in the call 
for rigour in this general area. In a typical scenario, we are provided with an 
infinity of random variables, indexed by the vertices of some graph such as the 
cubic lattice, and which have some dependence structure governed by the geometry 
of the graph. Thus mathematicians and physicists have had further cause to relate 
probability and geometry. One major outcome of the synthesis of ideas from physics 
and probability is the theory of Gibbs states, [59], which is now established as a 
significant branch of probability theory. 

A classic example of a Gibbs state is the (Lenz-)Ising model [89] for a ferromag- 
net. When formulated on the bounded region B of the square lattice L^, a random 
variable ax taking values —1 and is assigned to each vertex x of S, and the 
probability of the configuration a is proportional to exp(— /3i7(cr)), where /3 > 
and the 'energy' Hia) is the negative of the sum of ax(7y over all neighbouring 
pairs X, y of B. This 'starter model' has proved extraordinarily successful in gener- 
ating beautiful and relevant mathematics, and has been useful and provocative in 
the mathematical theory of phase transitions and cooperative phenomena (see, for 
example, [50]). 

There are many possible generalisations of the Ising model in which the (Jx may 
take a general number q of values, rather than g = 2 only. One such generalisation, 
the so-called Potts model [126], has attracted especial interest amongst physicists, 
and has displayed a complex and varied structure; for example, when q is large, 
it enjoys a discontinuous phase transition, in contrast to the continuous transition 
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believed to take place for small q. Ising/Potts models are the first of three princi- 
pal ingredients in the story of random-cluster models. Note that they arc 'vertex 
models' in the sense that they involve random variables ax indexed by the vertices 
X of the underlying graph. 

The '(bond) percolation model' was inspired by problems of physical type, and 
emerged from the mathematics literature of the 1950s [29, 150]. In this model for a 
porous medium, each edge of a graph is declared 'open' (to the passage of fluid) with 
probability p, and 'closed' otherwise, different edges having independent states. The 
problem is to determine the typical large-scale properties of connected components 
of open edges, as the parameter p varies. Percolation theory is now a mature part of 
probability, at the core of the study of random media and interacting systems, and 
it is the second ingredient in the story of random-cluster models. Note that bond 
percolation is an 'edge model', in that the random variables are indexed by the set 
of edges of the underlying graph. [There is a variant termed 'site percolation' in 
which the vertices are open/closed at random rather than the edges.] 

The third and final ingredient preceded the first two, and is the theory of elec- 
trical networks. Dating back at least to the 1847 paper [102] of Kirchhoff, this sets 
down a method for calculating macroscopic properties of an electrical network in 
terms of its local structure. In particular, it explains the relevance of counts of cer- 
tain types of spanning trees of the graph. In the modern vernacular, an electrical 
network on a graph G may be studied via the properties of a 'uniformly random 
spanning tree' on G (see [16]). 

These three ingredients seemed fairly distinct until Fortuin and Kasteleyn dis- 
covered, around 1970, [53, 54, 55, 94], that each features in a certain way within 
a family of probability measures of 'edge models', paramcterised by two quanti- 
ties, p G [0,1] and q G (0, cxo). [In actuality, electrical networks arise as a weak 
limit of such measures.] These models they termed 'random-cluster models', and 
they developed the basic theory — correlation inequalities and the like — in a 
series of papers published thereafter. The true power of random-cluster models as 
a mechanism for studying Ising/Potts models has emerged only gradually over the 
intervening thirty years. 

We note in passing that the genesis of the random-cluster model lay in Kaste- 
leyn's observation that each of the three ingredients above satisfies certain se- 
ries/parallel laws: any two edges in series (or parallel) may be replaced by a single 
edge in such a way that, if the interaction function is adapted accordingly, then the 
distributions of large-scale properties remain unchanged. 

The family of random-cluster measures (that is, probability measures which gov- 
ern random-cluster models) is not an extension of the Potts measures. The rela- 
tionship is more sophisticated, and is such that correlations for Potts models cor- 
respond to connections in random-cluster models. Thus the correlation structure 
of a Potts model may be studied via the stochastic geometry of a corresponding 
random-cluster model. The intuition behind this geometrical study comes often 
from percolation, of which the random-cluster model is indeed an extension. 

It turns out that, in many situations involving ferromagnetic Ising/Potts models, 
the best way forward is via the random-cluster model. As examples of this we 
mention the existence of discontinuous phase transitions [105], exact computations 
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in two dimensions [14] , the verification of the Wulff construction for Ising droplets 
[36], and the Dobrushin theory of interfaces and non-translation-invariant measures 
[62]. As a major exception to the mantra 'everything worth doing for Ising/Potts 
is done best via random-cluster', we remind the reader of the so-called random- 
current expansion for the Ising model, wielded with effect in [1, 2, 5] and elsewhere. 
The random-current method appears to be Ising-specific, and has enabled a deep 
analysis of the Ising model unparalleled in more general Potts systems. (See Section 
5.3.) 

The primary target of this review is to summarise and promote the theory of 
random-cluster models for its own sake. In doing so, we encounter many results 
having direct impact on Ising/Potts systems, but we shall not stress such connec- 
tions. Some of the theory has been discovered several times by apparently inde- 
pendent teams; whilst making a serious attempt to list key references, we apologise 
for unwitting omissions of which there will certainly be a few. The large number 
of references to work of the author is attributable in part to the fact that he is 
acquainted with these contributions. 

It is a lesson in humility to return to the original Fortuin-Kasteleyn papers [53, 
55], and especially [54], where so much of the basic theory was first presented. These 
authors may not have followed the slickest of routes, but they understood rather 
well the object of their study. Amongst the many papers of general significance 
since, we highlight: [4], which brought the topic back to the fore; [47], where the 
coupling between Potts and random-cluster models was so beautifully managed; 
[69], where the random-cluster model was studied systematically on infinite grids; 
[81], which links the theory to several other problems of interest in statistical me- 
chanics; and [61], where random-cluster models are placed in the perspective of 
stochastic geometry as a tool for studying phase transitions. 

This review is restricted mostly to core material for random-cluster models on 
the nearest-neighbour cubic lattice in a general number d of dimensions. Only in 
passing do we mention such subjects as extensions to long-range systems [4], mean- 
field behaviour in high dimensions [99], and mixing properties [7]. Neither do we 
stress the impact that graphical methods of the random-cluster type have had on a 
variety of other disordered systems, such as the Ashkin-Teller model [12, 81, 122, 
132, 148], the Widom-Rowlinson model [39, 40, 60, 81, 146], or on methods for 
simulating disordered physical systems [41, 42, 141, 149]. 

Random-cluster methods may be adapted to systems with random interactions 
[3, 72, 81], and even to non-ferromagnetic systems of Edwards- Anderson spin-glass 
type [46, 117, 118] where, for example, they have been used to prove that, for a given 
set { Je} of positive or negative interactions, uniqueness of the infinite-volume Gibbs 
measure for the ferromagnetic system having interactions { | </e | } implies uniqueness 
for the original system. 

Amongst earlier papers on random-cluster models, the following include a degree 
of review material: [4, 23, 61, 66, 81, 113]. 

Notwithstanding the fairly mature theory which has evolved, there remain many 
open problems including some of substantial significance for the general area. Many 
of these are marked in the text with the acronym OP. 
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2. Potts and random-cluster processes 



We write n{f) for the expectation of a random variable / under a probability 
measure /x. 

2.1 Random-cluster measures 

Let G = {V, E) he a finite graph. An edge e having end vertices x and y is written 
as e = {x,y). A random-cluster measure on G is a member of a certain class of 
probability measures on the set of subsets of the edge set E. We take as state space 
the set Q = {0, 1}"^, members of which are vectors lv = (w(e) : e E E). We speak 
of the edge e as being open (in iv) if a;(e) = 1, and as being closed if uj{e) = 0. For 
a; e n, let r]{u>) — {e & E : u>{e) = 1} denote the set of open edges, and let k{uj) 
be the number of connected components (or 'open clusters') of the graph (V, r]{uj)). 
Note that k{uj) includes a count of isolated vertices, that is, of vertices incident to 
no open edge. We assign to Q the u-field of all its subsets. 

A random- cluster measure on G has two parameters satisfying < p < 1 and 
g > 0, and is the measure (f)p^q on the measurable pair (fi, J^) given by 



This measure differs from product measure only through the inclusion of the term 
q'^i^) , Note the difference between the cases q < 1 and q > 1: the former favours 
fewer clusters, whereas the latter favours many clusters. When q = 1, edges are 
open/closed independently of one another. This very special case has been studied 
in detail under the titles percolation and random graphs; see [25, 71, 90]. Perhaps 
the most important values of q are the integers, since the random-cluster model 
with q e {2,3,...} corresponds, in a way sketched in the next two sections, to 
the Potts model with q local states. The bulk of this review is devoted to the 
theory of random-cluster measures when q >1. The case q < I seems to be harder 
mathematically and less important physically. There is some interest in the limit 
as g I 0; see Sections 2.4 and 3.6. 

We shall sometimes write 4>G,p,q for 4'p,q^ when the choice of graph G is to be 
stressed. Samples from random-cluster measures on 1? are presented in Figure 2.1. 

2.2 IsiNG and Potts models 

In a famous experiment, a piece of iron is exposed to a magnetic field. The field 
is increased from zero to a maximum, and then diminishes to zero. If the temper- 
ature is sufficiently low, the iron retains some 'residual magnetisation', otherwise 
it does not. There is a critical temperature for this phenomenon, often called the 
Curie point. The famous (Lenz-)Ising model for such ferromagnetism, [89], may be 
summarised as follows. One supposes that particles are positioned at the points of 




where the 'partition function', or 'normalising constant', Z is given by 
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p = 0.51 p = 0.55 p = 0.70 



Figure 2.1. Samples from the random-cluster measure with q = 1 on the box [0, 40]^ 
of the square lattice. We have set q = 1 for ease of programming, the measure being 
of product form in this case. The critical value is Pc(l) = \- Samples with more 
general values of q may be obtained by the method of 'coupling from the past', as 
described in Section 8.2. 

some lattice in Euclidean space. Each particle may be in either of two states, rep- 
resenting the physical states of 'spin up' and 'spin down'. Spin-values are chosen at 
random according to a certain probability measure, known as a Gibbs state, which 
is governed by interactions between neighbouring particles. This measure may be 
described as follows. 

Let G = {V, E) be a finite graph. We think of each vertex v e F as being 
occupied by a particle having a random spin. Since spins are assumed to come in 
two basic types, we take as sample space the set S = {—1, +1}^. The appropriate 
probability mass function X/3,j,h on S has three parameters satisfying < /3, J < oo 
and /i G M, and is given by 

A;3,j,/.(a) = ^e-^^(-), aeE, 

where the partition function Zi and the 'Hamiltonian' : S — > R are given by 

o-GE e={x,y)eE xEV 

The physical interpretation of P is as the reciprocal 1/T of temperature, of J as 
the strength of interaction between neighbours, and of h as the external field. For 
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reasons of simplicity, we shall consider only the case of zero external field, and we 
assume henceforth that h = 0. Each edge has equal interaction strength J in the 
above formulation. Since f3 and J occur only as a product /3J, the measure A^^j^o 
has effectively only a single parameter. In a more complicated measure not studied 
here, different edges e are permitted to have different interaction strengths Jg. 

As pointed out by Baxter, [14], the Ising model permits an infinity of general- 
isations. Of these, the extension to so-called 'Potts models' has proved especially 
fruitful. Whereas the Ising model permits only two possible spin-values at each 
vertex, the Potts model [126] permits a general number q E {2,3, . . .}, and is given 
as follows. 

Let q be an integer satisfying g > 2, and take as sample space S = {l,2,...,g}^. 
Thus each vertex of G may be in any of q states. The relevant probability measure 
is now given by 

where Zp is the appropriate normalising constant, 

and 6u,v is the Kronecker delta. When q = 2, we have that Sa^,(Ty = |(1 + c^Cj/)- 
It is now easy to see in this case that the ensuing Potts model is simply the Ising 
model with an adjusted value of J. 

2.3 Random-cluster and Ising-Potts coupled 

It was Fortuin and Kasteleyn [53, 54, 55, 94] who discovered that Potts models 
may be recast as random-cluster models, and furthermore that the relationship 
between the two systems facilitates an extended study of phase transitions in Potts 
models. Their methods were essentially combinatorial. In the more modern game, 
we construct the two systems on a common probability space, and then observe 
their relationship through their realisations. There may in principle be many ways 
to do this, but the standard coupling reported in [47] is of special value. 

Let g G {2, 3, . . . }, < p < 1, and let G = {V, E) be a finite graph, as before. We 
consider the product sample space E x O where E = {1, 2, . . . , q}^ and VL = {0, 1}^ 
as above. We now define a probability mass function on E x Q by 

ll{a,u) OC Jj|(l -p)(^a;(e),0 +P<^u;(e),l<^e((^)}, (o", ^) G E X J], 

e&E 

where (5e(cr) = Scj^^ay for e = {x,y) e E. Elementary calculations reveal the 
following facts. 

(a) Marginal on E. The marginal measure jJ-iia) = X^i^^n is given by 

lii{a) (X ex.pl (3 ^ J5e(o-) > 

^ e ^ 

where p=l — e~^^ . This is the Potts measure. 
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(b) Marginal on Q. The second marginal of /i is 



o-GE ^ e 



This is the random-cluster measure. 

(c) The conditional measures. Given cu, the conditional measure on E is obtained 
by putting (uniformly) random spins on entire clusters of u (of which there 
are k{uj)). These spins are constant on given clusters, and are independent 
between clusters. Given cr, the conditional measure on Q is obtained by 
setting uj{e) = if 5e(cr) = 0, and otherwise a;(e) = 1 with probability p 
(independently of other edges). 

In conclusion, the measure |U is a coupling of a Potts measure i^p^j^q on V, 
together with the random-cluster measure 4>p^q on Vt. The parameters of these 
measures are related by the equation p=l — e~^'^ . Since < p <1, this is possible 
only if pj > 0. 

This special coupling may be used in a particularly simple way to show that cor- 
relations in Potts models correspond to open connections in random-cluster models. 
When extended to infinite graphs, this implies as discussed in Section 4 that the 
phase transition of a Potts model corresponds to the creation of an infinite open 
cluster in the random-cluster model. Thus arguments of stochastic geometry, and 
particularly those developed for the percolation model, may be harnessed directly 
in order to understand the correlation structure of the Potts system. The basic 
step is as follows. 

We write {x ^ y} for the set of all a; G O for which there exists an open path 
joining vertex x to vertex y. The complement of the event {x ^ y} is denoted 
{x ^ y}. 

The 'two-point correlation function' of the Potts measure 77/3 on the finite 
graph G = (y, E) is defined to be the function rp^j^q given by 



The term q is the probability that two independent and uniformly distributed 
spins are equal. The 'two-point connectivity function' of the random-cluster mea- 
sure (f)p^q is defined as the function (t>p,q{x ^ y) for x,y eV, that is, the probability 
that X and y are joined by a path of open edges. It turns out that these 'two-point 
functions' are (except for a constant factor) the same. 

Theorem 2.1 (Correlation/connection) [94]. If q e {2,3,...} and p = 1 — 

e~^'^ satisfies < p < 1, then 



rf3,J,q{x, y) = T^I3,J,q{(7x = CTy) 



1 



x,y eV. 



Tf3,J,q{x,y) ^ {1- q ^)(j)p,q{x^ 



y), 



x,y ^ V. 
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Proof. The indicator function of an event A is denoted 1a- We have that 

ri3,J,qix, y) = ^[^{a^=ay}{<y) " ^) 

<T,U> 

LO a 

= (1 - Q.'^)^p,q{^ ^ y), 

where /x is the above couphng of the Potts and random-cluster measures. □ 

The theorem may be generahsed as foUows. Suppose wc are studying the Potts 
model, and are interested in some 'observable' / : S ^ E. The mean value of f{a) 
satisfies 

T^(3,J,q{f) = ^f{(^)T^I3,J,q{(7) = ^ f {(t) n{a , u) 
(J a,(jj 

= Y,Fi^)^P,qi^) = ^P,qiF) 

where F : f2 ^ R is given by 

F(a;)=M/l^) = E/HM^I^)- 

The above theorem is obtained in the case /(cr) = 5a^,cTy — where x,y E V. 

The Potts models considered above have zero external field. Some complications 
arise when an external field is added; see the discussions in [8, 23]. 

2.4 The limit as J, 

Let G = (y, E) be a finite connected graph, and let (j)p^q be the random-cluster 
measure on the associated sample space O = {0, 1}^. We consider first the weak 
limit of (f)p^q as g J, for fixed p e (0, 1). This limit may be ascertained by observing 
that the dominant terms in the partition function 

Z{p,q) = 

wen 

are those for which A; (a;) is a minimum, that is, those with A; (a;) = 1. It follows that 
lim^j^o '^p,q is precisely the product measure (f)p^i (that is, percolation with intensity 
p) conditioned on the resulting graph {V, r]{u})) being connected. A more interesting 
limit arises if we allow p to converge to with q, as follows. 

The random-cluster model originated in a systematic study by Fortuin and 
Kasteleyn of systems of a certain type which satisfy certain parallel and series laws. 
Electrical networks are the best known such systems — two parallel (respectively, 
series) connections of resistances ri and r2 may be replaced by a single connection 
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with resistance (rf ""^ + r2^)~^ (respectively, ri + r2). Fortuin and Kasteleyn [55] 
realised that the electrical-network theory of a graph G is related to the limit as 
g I of the random-cluster model on G. Their argument may be expanded as 
follows. 

Suppose p = pq is related to q in such a way that p — > and q/p ^ as q ^ 0. 
We may write Z{p, q) as 

Note that — p) — ^ and ^(1 — — > as ^ 0. Now k{u)) > 1 and 
|?7(u;)| +k{u>) > \V\ for all u> E these two inequalities are satisfied simultaneously 
with equality if and only if 77(0;) is a spanning tree of G. It follows that, in the limit 
as g — > 0, the 'mass' is concentrated on such configurations, and it is easily seen 
that the limit mass is uniformly distributed. That is, lim^jo 4'p,q is a probability 
measure which selects, uniformly at random, a spanning tree of G; in other words, 
the limit measure is 0i conditioned on the resulting graph being a spanning tree. 

The link to the theory of electrical networks is now provided by Kirchhoff's 
theorem [102], which expresses effective resistances in terms of counts of spanning 
trees. See also [79]. 

The theory of random spanning trees is beautiful in its own right (see [16]), and 
is linked in an important way to the emerging field of stochastic growth processes 
of 'stochastic Lowner evolution' (SLE) type (see [111, 130]), to which we return in 
Section 6.4. Another limit emerges if p = q and q i 0, namely uniform measure on 
the set of forests of G. More generally, take p = aq where a e (0, 00) is constant, 
and take the limit as g | 0. The limit measure is the percolation measure (f)f3^i 
conditioned on the non-existence of open circuits, where (3 — a/{l + a). Ifp/g— >0 
as p, g — >^ 0, the limit measure is concentrated on the empty set of edges. 

2.5 Rank-generating functions 

The partition functions of Potts and random-cluster measures are particular evalua- 
tions of rank-generating functions, defined as follows. The rank-generating function 
of the simple graph G = (V, E) is the function 

Wg{u,v) = Yl u'^^^"^v^^^'\ u.veR, 

E'CE 

where r(G") = \V\ - k{G') is the rank of the graph G' = {V,E'), and c{G') = 
\E'\ — \V\ + k{G') is its co-rank-^ here, k{G') denotes the number of components 
of the graph G' . The rank-generating function has various useful properties, and 
occurs in several contexts in graph theory; see [20, 142]. It crops up in other forms 
also. For example, the function 

Tg{u, v) = {u- 1)\''\-^Wg{{u - v-l) 

is known as the dichromatic (or Tutte) polynomial, [142]. The partition function 
Z — Zq of the random-cluster measure on G with parameters p, q is easily seen to 
satisfy 

Z„^,l-l(l-p)l-IW.(^,^). 
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a relationship which provides a hnk with other classical quantities associated with 
a graph. See [20, 21, 53, 144] also. 

3. Infinite-volume random-cluster measures 

It is in the infinite-volume limit that random-cluster measures exhibit phase tran- 
sitions. There are two ways of constructing random-cluster measures on infinite 
graphs, namely by taking weak limits as a finite domain approaches the infinite 
system, and by studying measures on the infinite graph having the 'correct' con- 
ditional versions. Such matters are discussed in this section, which begins with a 
summary of certain valuable properties of random-cluster measures on finite graphs. 

3.1 Stochastic ordering 

The stochastic ordering of measures provides a technique fundamental to the study 
of random-cluster measures. Let G = {V, E) be a finite or countably infinite graph 
as above; let O = {0, l}"^, and let be the cr-field of O generated by the finite- 
dimensional cylinders. Note first that is a partially ordered set with partial order 
uji < u>2 if <^i(e) < (e) for all e. A random variable / : f2 — > R is called increasing 
if f{uJi) < /(<^2) whenever uji < u>2- An event A & T is called increasing if its 
indicator function 1a is increasing. The word 'decreasing' should be interpreted in 
the natural way. Given two probability measures fii, ^2 on Q, we write fii <st A*2 5 
and say that /xi is stochastically smaller than /X2, if A*i(/) < A*2(/) for all bounded 
increasing random variables f on Q. 

We return now to the case when G is a finite graph. Let ni , 112 be probability 
measures on Q, and assume for the moment that the fii are strictly positive in 
the sense that Hiito) > for all u; G O. An important sufficient condition for the 
inequality ni <st A*2 was found by HoUey [88], namely that 

//i(a;i Va;2)At2(wi A^2) > Ati(^i)At2(w2) for aU uji,uj2 e 9., 

where uji V UJ2 and lji A 0^2 are the maximum and minimum configurations given 
respectively as max{c<;i(e), 0^2(6)} and nim{uji{e) , uj2{e)} , for e & E. A probability 
measure on is said to have the FKG lattice property if 

IJ,{uJl y UJ2)lJ-{iO\ /\UJ2) > IJ,{u!l)fi{u!2) for all UJl,U!2 £ ^, 

amd it is a consequence of HoUey's argument that any strictly positive measure with 
the FKG lattice property satisfies the so-called FKG inequality. This amounts to 
the following for random-cluster measures. 

Theorem 3.1 (FKG inequality) [54, 56]. Suppose that < p < 1 and q > 1. 

If f and g are increasing functions on Q, then 4>p,q{fg) > 4)p^q{f)4)p^q{g). 

Specialising to indicator functions, we obtain that 

(f>p,q{A f] B) > (j)p^q{A)(j)p^q{B) for increasing events A^ B, 

whenever q>l. It is not difficult to see that the FKG inequality does not generally 
hold when < q < 1. 

HoUey's theorem leads easily to the following comparison inequalities, which 
were first proved by Fortuin. 
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Theorem 3.2 (Comparison inequalities) [54]. It is the case that 

(j)p',q' < (j)p,q if > g, q' > l, andp' < p, 
(f>p',g' > (f>p,g if q' >q, q'> 1, and > 

3.2 A DIFFERENTIAL FORMULA 

One way of estimating the probability of an event A is via an estimate of its 
derivative d(f)p^q{A)/dp. When q = 1, there is a formula for this derivative which 
has proved very useful, and which is commonly attributed to Russo, see [13, 71, 
131]. This formula may be generalised to random-cluster measures as follows. The 
proof is an exercise in the differentiation of summations. 

Theorem 3.3 [19]. Let < p < 1, q > 0, and let (pp^q he the corresponding 
random- cluster measure on a finite graph G = {V, E). Then 

^(Pp,q{A) = ^^^^{(I>pM'^a) - (f>p,q{\v\)(f>p,q{A)} 

for any event A, where \r]\ = \r]{u!)\ = X^ees'^l^) number of open edges of 

the configuration u. 



3.3 Conditional probabilities 

Whether or not an edge e is open depends on the configuration on i? \ {e}, and a 
further important property of random-cluster measures summarises the nature of 
this dependence. 

For e e E, we denote by G \ e (respectively, G.e) the graph obtained from G 
by deleting (respectively, contracting) e. We write fig = {0, l}-^^^^^; for a; e f2 we 
define e Qe by uJe{f) = <^(/) for f ^ e. For e = {x, y), we write Kg for the event 
that X and y are joined by an open path not using e. 



Theorem 3.4 [54]. Let e e E. We have that 
(l>G,p,q{uj\uj{e) = j) = 

and 



(t'G\e,p,q{^e) if j = 0, 
4>G.e,p,q{f^e) if 3 = 



p if(^e& Ke, 

^G,p,qHe) = 1 I a;e) = <! P 

p+ {1 - p)q 



ifUe i Ke. 



That is to say, the effect of conditioning on the absence or presence of an edge 
e is to replace the measure 4>G,p,q by the random-cluster measure on the respective 
graph G \ e or G.e. Secondly, the conditional probability that e is open, given the 
configuration elsewhere, depends only on whether or not K^. occurs, and is then 
given by the stated formula. The proof is elementary. The final equation of the 
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theorem leads to properties of random-cluster measures referred to elsewhere as 
'insertion tolerance' and the 'finite-energy property'. 

3.4 Infinite- VOLUME weak limits 

In studying random-cluster measures on infinite graphs, we restrict ourselves to the 
case of the hypercubic lattice in d dimensions, where d > 2; similar observations 
are valid in greater generality. Let d>2, and let 1/^ be the set of all d-vectors of 
integers; for x e Z*^, we normally write x = {xi,X2, ■ ■ ■ , Xd). For x,y & Z*^, let 



i=l 



We place an edge {x, y) between x and y if and only if ||a; — y|| = 1; the set of such 
edges is denoted by E*^, and we write = (Z*^, E*^) for the ensuing lattice. For any 
subset S oiZ^, we write dS for its boundary, that is, 

dS = {seS: {s, t) e E^ for some t ^ S}. 

Let Q = {0, 1}'^'*, and let be the u-field of subsets of Q generated by the finite- 
dimensional cylinders. The letter A is used to denote a finite box of Z^, which 
is to say that A = IliLii^j'^/i] some x,y G Z'^; we interpret [xi^yi] as the set 
{xi, Xi + l,Xi + 2, . . . ,yi}. The set A generates a subgraph of L*^ having vertex set 
A and edge set Ea containing all (x, y) with x, y e A. 

We are interested in the 'thermodynamic limit' (as A f Z'^) of the random-cluster 
measure on the finite box A. In order to describe such weak limits, we shall need 
to introduce the notion of a 'boundary condition'. 

For ^ G O, we write for the (finite) subset of O containing all configurations 
uj satisfying uj{e) = ^(e) for e G E'^ \ Ea; these are the configurations which 'agree 
with ^ off A'. For ^ G Q and values of p, q satisfying 0<p<l,(j>0, we 
define to be the random-cluster measure on the finite graph (A, Ea) 'with 

boundary condition ^'; this is the equivalent of a 'specification' for Gibbs states. 
More precisely, let 0a p g ^® probability measure on the pair (il, T) given by 

-^\\[p-^^\l-pY--^Aq^^-^^) ifcG^, 

otherwise, 

where k{u!, A) is the number of components of the graph (Z'^, 'r]{uj)) which intersect 
A, and where Z^p q is the appropriate normalising constant 

4.,.= E { ^^''^^''^(l-^')'"^^^4^'^"^^• 
Note that 4>i,p,q{^i) = 1- 
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Definition 3.5. Let < p < 1 and q > 0. A probability measure </> on (fi,^) is 
called a limit random- cluster measure with parameters p and g if there exist ^ G O 
and a sequence A = (A^ : n > 1) of boxes satisfying A„ ^ Z*^ as n — > oo such that 

A^,P,q as n ^ OO. 

The set of all such measures (f) is denoted by yVp,q, and the closed convex hull of 
yVp,q is denoted coWp,g. 

In writing A^ — > Z'^ we mean that, for all m, A^ ^ [— m, m]^ for all large n. The 
arrow '=^' denotes weak convergence. 

It might seem reasonable to define a limit random-cluster measure to be any 
weak limit of the form lim^^oo (^a^^ p q for some sequence (^n : n > 1) of members 
of Q. and some sequence A = (A„ : n > 1) of boxes satisfying A^ ^ U^. It may 
however be shown that this adds no extra generality to the class as defined above, 
[69]. The dependence of the limit measure (f) on the choice of sequence (A^) can be 
subtle, especially when q < 1 (OP). 

It is standard that yVp,q 7^ for all < p < 1, g > 0, and one way of seeing this 
is as follows. The sample space Q, is the product of discrete spaces, and is therefore 
compact. It follows that any class of probability measures on O is tight, and hence 
relatively compact (see the account of Prohorov's theorem in [22]), which is to say 
that any infinite sequence of probability measures contains a weakly convergent 
subsequence. 

When does the limit lim^^oo </'a„ p q exist, and when does it depend on the choice 
of boundary condition ^? The FKG inequality provides a route to a partial answer 
to this important question. Suppose for the moment that q > 1. Two extremal 
boundary conditions of special importance are provided by the configurations 
and 1, comprising 'all edges closed' and 'all edges open' respectively. One speaks 
of configurations in as having 'free' boundary conditions, and configurations in 
as having 'wired' boundary conditions. 

Theorem 3.6 (Thermodynamic Umit) [4, 27, 54, 66, 69]. Suppose < p < 1 

and q>l. 

(a) Let A = (A^ : n > 1) be a sequence of boxes satisfying A„ ^ Z'' as n — > 00. 
The weak limits 

4>i,q = }}^Ar,,p,q, for b = 0,1, 

exist and are independent of the choice of A. 

(b) We have that each (j)^ ^ is translation-invariant, and 

(j)l,q <st (P <st (l>l,q for all (f) e Wp^q. 

(c) For b = 0,1, the measure (j)^ ^ is ergodic, in that any translation-invariant 
random variable is (j)^ q-a.s. constant. 

The FKG inequality underlies all parts of Theorem 3.6. The claim (c) of ergod- 
icity has until recently been considered slightly subtle (see the discussion after the 
forthcoming Theorem 3.9) but an easy proof may be found in [113]. 
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It follows from the inequality of part (b) that | Wp,g| = 1 if and only if 4)^ ^ = 0^ g. 
It is an important open problem to determine for which p, q this holds, and we shall 
return to this question in Section 5 (OP). For the moment, we note one sufficient 
condition for uniqueness, proved using a certain convexity property of the logarithm 
of a partition function Z. 

Theorem 3.7 [67, 69]. Let q > I. There exists a subset Vg of [0,1], at most 
countably infinite in size, such that (p^^^ = 4>l^q, and hence \yVp,q\ = 1, if P ^T^q- 

It is believed but not proved (OP) that: for any given Q > 1, "Dq either is empty 
or consists of a singleton (the critical point, to be defined in Section 4), the former 
occurring if and only if q is sufficiently small. 

3.5 Random-cluster measures on infinite graphs 

One may define a class of measures on the infinite lattice without having recourse 
to weak limits. The following definition of a random-cluster measure is based upon 
the Dobrushin-Lanford-Ruelle (DLR) definition of a Gibbs state, [44, 59, 106]. It 
was introduced in [66, 67], and discussed further in [27, 69]. For any box A, we 
write 7a for the cr-field generated by the set {a;(e) : e e E'^ \ Ea} of states of edges 
having at least one end vertex outside A. 

Definition 3.8. Let < p < 1 and g > 0. A probability measure on (fi,.F) is 
called a random- cluster measure with parameters p and q if 

for all ^ e and all finite boxes A, (/)(^ | 7^)(^) = (/»^ p^g(A) for 0-a.e. ^. 

The set of such measures is denoted Ttp^q- 

The condition of this definition amounts to the following. Suppose we are given 
that the configuration off the finite box A is that of ^. Then, for almost every 
^ e the (conditional) measure on A is simply the random-cluster measure with 
boundary condition ^. No further generality is gained by replacing the finite box A 
by a general finite subset of U^. 

Some information about the structure of 7^p,g, and its relationship to Wp.g, is 
provided in [69]. For example, for all p, g, 7^p,q is non-empty and convex. We have 
no proof that Wp,g C 7^p,q, but we state one theorem in this direction. For a; G fi, 
let I{uj) be the number of infinite open clusters of uj. We say that a probability 
measure on (fi, .F) has the Q /I- infinite- cluster property if 0(7 e {0, 1}) = 1. 

Theorem 3.9 [69, 70, 73]. Let Q < p < I and q > If (f) e co Wp,q and <j) has 
the Q/l-infinite- cluster property, then (j) £ T^p,q- 

Since, [30], any translation- invariant probability measure satisfying a finite- 
energy property (see the discussion after Theorem 3.4) necessarily has the 0/1- 
infinite-cluster property, we have that all translation-invariant members of co Wp,g 
lie in TZp^q. Suppose for the moment that q > 1. By Theorem 3.6(b), the weak 
limits 6 = 0, 1, are translation-invariant, and therefore they belong to 7^p,q. It 
is not difficult to see, by the FKG inequality, that 

(*) 4>l,q <st <st for all (j) e 7^p,g, 
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and it follows that \'R^p,q\ = 1 if and only if 0°^^ = 0^ g. The claim of ergodicity 
in Theorem 3.6(c) is one consequence of the extremality (*) of the 0^ ^ within the 
class Tlp^q (see also [113, page 1113]). 

It may be seen by an averaging argument, [69], that coWp^g necessarily con- 
tains at least one translation-invariant measure, for all p e [0, 1] and q e (0, oo). 
Therefore, 7^p,g is non-empty for all p and q. 

We note that Theorem 3.9, and particularly the 0/1-infinite-cluster property, is 
linked to the property of so-called 'almost sure quasilocality', a matter discussed in 
[121]. 

3.6 The case q<l 

The FKG inequality, a keystone of many arguments when g > 1, is not valid 
when q < 1. Consequently, many fundamental questions are unanswered to date, 
and the theory of random-cluster models on a finite graph G — {V, E) remains 
obscure when g < 1. The intuition is that certain positive correlations should be 
replaced by negative correlations; however, the theory of negative correlation is 
more problematic than that of positive correlation (see [120]). We return to this 
point later in this subsection. 

As referred to above, there is an existence proof of infinite-volume weak limits 
and random-cluster measures for all g > 0. On the other hand, no constructive proof 
is known of the existence of such measures when q < 1 (OP). More specifically, 
the existence of the weak limits limy^^^d 4>\p q, b — 0, 1, is not known when q < 1. 
The best that can be shown currently is that the two limits exist and are equal 
when p is either sufficiently small or sufficiently large, [73] . This may be achieved 
by comparison with percolation models having different values of p, very much as 
in [69] (the claim for small p may also be shown by the arguments of [49, 51]). 

The theory of percolation gives a clue to a possible way forward. When q = 1, 
the FKG inequality is complemented by the so-called 'disjoint-occurrence' (or 'BK') 
inequality. This latter inequality is said to be valid for a measure /i if fi{A o 
B) < n{A)n{B) for all increasing events A, B, where A o B is the event that A 
and B occur disjointly (see [18, 71] for a discussion of this and the more general 
'Reimer inequality' [129]). The disjoint-occurrence inequality has been established 
for classes of measures which are only slightly more general than product measures, 
and it is an interesting open question whether it is valid for a wider class of measures 
of importance (OP). It has been asked whether the disjoint-occurrence inequality 
could be valid for random-cluster measures with g < 1 (OP). A positive answer 
would aid progress substantially towards an understanding of limit random-cluster 
measures. 

We illustrate this discussion about disjoint-occurrence with the following test 
question (OP): is it generally the case that the random-cluster measure 0p^q on G 
satisfies 

(*) (edges e and / are open) < 4>p,q{e is open)(/)p^g(/ is open) 

ior f and g < 1? (See [120].) This equation would be a very special instance of 
the disjoint-occurrence inequality. A further restriction arises if we take the limit 
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as J, 0; recall the discussion of Section 2.4. This leads to certain open questions of 
a purely graph-theoretic type, which combinatorial theorists might elevate to the 
status of conjectures. The first such question is the following. Let K(ei, 62, . . . ) 
be the number of subsets F of the edge set E, containing ei, 62, . . . , such that the 
graph (y, F) is connected. Is it the case that (OP) 

(**) K{e,f)K{0)<K{e)K{f) if e ^ /? 

(See [93].) In the second such question, we ask if the same inequality is valid with 
K{ei, 62, ... ) redefined as the number of subsets F containing ei, 62, . . . such that 
{V,F) is a forest (OP). These two questions are dual to one another in the sense 
that the first holds for a planar graph G if and only if the second holds for its planar 
dual. Explicit computations have confirmed the forest conjecture for all graphs G 
having nine or fewer vertices, [78]. 

In the 'intermediate regime', with K{ei, 62, ... ) redefined as the number of span- 
ning trees (that is, connected forests) of G containing ei, 62, . . . , the corresponding 
inequality is indeed valid. An extra ingredient in this case is the link to electri- 
cal networks, and particularly the variational principle known as the Thomson or 
Dirichlet principle (see [45]). Further results and references are provided in [16]. 
Substantially more is known for spanning trees, namely a general result concerning 
the 'negative association' of the uniform measure on the set of spanning trees of G, 
[48]. 

We note a more general version of conjecture (**), namely 

K^{e, f)K^{0) < K^{e)KM) for e ^ /, < a < 00, 

where 

Ka(ei,e2, . . .) = ^ a'^l. 

FCE 
FD{ei,e2,...} 
{V,F) connected 

This is equivalent to (*) in the limit as g | 0, where a = p/(l — p). 

By other means one may establish a certain non-trivial monotonicity when q < 1, 
but by a more complicated reasoning than before involving a property of convexity 
of the logarithm of the partition function. Namely, the mean number of open edges 
is non-decreasing in p, for < g < 00, [69]. 

4. Phase transition, the big picture 

Phase transition in a Potts model corresponds to the creation of an infinite open 
cluster in the corresponding random-cluster model. There are rich predictions con- 
cerning the nature of such a phase transition, but these have been proved only in 
part. This section is a summary of the expected properties of the phase diagram for 
different dimensions d and cluster-weighting factors q. The corresponding rigorous 
theory is described in Sections 5 and 6. 

4.1 Infinite open clusters 
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We assume henceforth that g > 1, and we concentrate here on the extremal random- 
cluster measures (f)^ ^ and (p^ ^. The phase transition of a random-cluster measure 
is marked by the onset of an infinite open cluster. We write {0 ^ oo} for the event 
that the origin is the endvertex of some infinite open path, and we define the 0^^^ 
percolation probability by 

e\p,q) = (f>l^{0^oo), 6 = 0,1. 

It is almost immediate by a stochastic-ordering argument that 9^{p,q) is non- 
decreasing in p, and therefore 

t > iip>pl{q), 
for critical points p^{q) given by 

p,^(g)=sup{p:^\p,g) = 0}, 6 = 0,1. 
It is an easy exercise to show that the number / of infinite open clusters satisfies: 



p,<i 



(/>1) 



if^^(p,g) = 0, 

1 if e^{p, q) > 0. 



We shall see in Section 5.2 that any infinite open cluster is 0^ ^-a.s. unique whenever 
it exists. 

We have by Theorem 3.7 that ^ = (p^ ^ for almost every p, whence q) = 
9^{p,q) for almost every p, and therefore p^iq) = pliq)- Henceforth we use the 
abbreviated notation Pc{q) = Pc{q) = pliQ)i we refer to Pc{q) as the critical 
point of the corresponding random-cluster measures. The non-triviality of Pc{q) 
may be proved by comparisons of random-cluster measures with product measures 
via Theorem 3.2. Recall the fact, [71, Chapter 1], that < pc{l) < 1 if d > 2. 

Theorem 4.1 [4]. We have for q >1 that 

gPc(l) 



Pc{l) <Pc{q) < 



l + (g-l)pc(l) 



When q is an integer satisfying q >2, the phase transition of the random-cluster 
model corresponds in a special way to that of the Potts model with the same value 
of q. An indicator of phase transition in the Potts model is the 'magnetisation', 
defined as follows. Consider a Potts measure n\ on A having parameters /3, J, g, 
and with '1' boundary conditions, which is to say that all vertices on the boundary 
dA are constrained to have spin value 1. Let ta = '^Xio'o = ^) — q~^, a quantity 
which represents the net effect of this boundary condition on the spin at the origin. 
The corresponding random-cluster measure (f)\ has parameters p = 1 — e~^'^ and q, 
and has wired boundary condition. We apply Theorem 2.1 to the graph obtained 
from A by identifying all vertices in dA, and we find that 

TA = (1 - g-')0i(O ^ dA). 
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The hmit function r = hmyv-|-z<* Ta. is called the magnetisation, it is a non-decreasing 
function of /3J and satisfies 

{=0 if (3 J is small, 
> otherwise. 

It is not hard to show, [4], that (/>\(0 ^ dK) ^ oo) as A | Z'^, whence 

T = (1 — ?~^)^^(P5 ?) wherep = l — e~^'^. Therefore there is long-range order in the 
Potts model (that is, r > 0) if and only if the origin lies in an infinite open cluster 
with strictly positive (/>p ^-probability In particular, Pc{q) = 1 — e~^'''^ where (3^ is 
the critical value of (3 for the Potts model in question. 

4.2 First- and second-order phase transition 

There is a rich physical theory of phase transitions in percolation, Ising, and Potts 
models, some of which has been made rigorous in the context of the random-cluster 
model. There follows a broad sketch of the big picture, a full rigorous verification of 
which is far from complete. Rigorous mathematical progress is described in Section 
5. 

I. The subcritical phase, p < Pc{q)- 
It is standard, [4], that 

cf>l,, = cl>l, iie'{p,q) = 0, 

implying that there exists a unique random-cluster measure whenever 9^{p, q) = 0. 
In particular, |Wp,g| = \T^p,q\ — 1 when < p < Pc{q)- Assume for the moment 
that p < Pc{q)i and denote the unique random-cluster measure by (t)p,q- By the 
definition of the critical point, all open clusters s. finite. It is believed 

that they have a tail which decays exponentially, in that there exist = ^{p, q) > 
and r] = r]{p,q) > such that 

0p,,(|C| = n) = e-Tn(i+o(i)), </,p^g(rad(C) = n) = e-^"(i+o(i)), as n ^ oo, 

where C denotes the open cluster containing the origin, and its radius rad(C) is 
defined as sup{||x|| : x e C}. Such exponential decay would be the starting point for 
a complete exploration of the subcritical phase. More detailed asymptotics should 
then emerge, including the Ornstein-Zernike decay of the connectivity functions: 

0p,g(O ^ en) ~ ^^^-""'^^"'''^ as n ^ OO, 

where = (n, 0, 0, . . . , 0) and ^{p, q) is termed the 'correlation length'. 

II. The supercritical phase, p > Pc{q)- 

This phase is characterised by the existence of one or more infinite open clusters 
(exactly one, in fact, for translation-invariant measures at least, see Section 5.2). 
It is believed that, as in the subcritical phase, we have that (pp ^ = ^ when 
p > Pciq)', this remains unproven in general. Thus the first main problem is to 
prove that there is a unique random-cluster measure when p > Pc (q) ■ 
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The theory of percolation, [71] , suggests a route towards understanding the ge- 
ometry of the supercritical phase, namely by developing a rigorous block renor- 
malisation argument. This should permit the use of theory developed when p is 
close to 1 in order to understand the model when p is close to Pc{q)- In particular, 
one expects an exponential estimate for the decay of the probabilities of long-range 
connections within finite open clusters, and a Wulff construction for the shape of 
such clusters. 

III. Near the critical point, p ~ Pcig)- 

The main open problem is to understand the way in which the nature of the phase 
transition depends on the value of q. It is believed that the transition is continuous 
and governed by critical exponents and scaling theory when q is small, and is 
discontinuous when q is large. Presumably there exists a threshold for q which 
separates the so-called 'second-order' (or continuous) transition from the so-called 
'first-order' (or discontinuous) transition. More specifically, it is believed that there 
exists Q = Q{d) satisfying 



(i) Assume that q < Q. 

- For any p, there exists a unique random-cluster measure, denoted (f)p,q- In 



- 9{p, q) = (j)p,q{0 ^ oo) is a continuous function of p. There is no percolation 

at the critical point, in the sense that 0{pc{q), q) = 0. 

- The edge-density h{p, q) — 4>p,q{e is open), viewed as a function of p, is con- 

tinuous at the critical point p — Pc{q)- [The letter e denotes a typical edge 
of the lattice.] 

- These functions and others have power-law singularities at Pc{q), and the 

associated critical exponents satisfy the scaling relations (see [71, Chapter 



- When d is large (how large depends on the value of q), these critical exponents 

take on their 'mean-field' values, and depend no further on the value of d. 

- There is no 'mass gap', in the sense that the correlation length ^(p, g) satisfies 

linipTp,((,)^(p, q) = oo. 

- Universality reigns, in that the critical exponents depend on the number d of 

dimensions but not on the choice of lattice. For example, the exponents 
associated with the square lattice are expected to be the same as those for 
the triangular lattice. 

- Assume d — 2 and 1 < g < 4. The process with p = Pc{q) converges as 

the lattice spacing shrinks to zero, the limit process when suitably defined 
being a stochastic Lowner evolution SLE^ having parameter k satisfying 
cos(47r/«;) = -^^/q, k e (4,8) (see Section 6.4 and [130]). 
(ii) Assume that q > Q. 




such that the following hold. 



particular 



Pc{q),q 



Pc{q),q 



9]). 
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- There exists a unique random-cluster measure if and only if p 7^ Pcig)- When 

d = 2 and p = Pciq), there are exactly two extremal members of 7^p,g, 
namely the free and the wired measures (pp q, b = 0, 1. When d > 3 and 
P — Pc{q) there exist other extremal members of TZp^q including a variety 
of non-translation-invariant measures. 

- We have that 0^{pc{q),q) = but 0^{pc{q),q) > 0. 

- The edge-density h{p^ q) is a discontinuous function of p at the critical point 

Pc{q)- 

- There is a 'mass gap' in the sense that the correlation length ^(p, q) satisfies 

limprpc(Q) 9) < 00. 

5. General results in d (> 2) dimensions 

The properties of the random-cluster model depend pivotally on whether the process 
is subcritical {p < Pc{q)), supercritical {p > Pc{q)), or critical {p ~ Pciq))- We 
consider these situations in turn, in each case identifying major results and open 
problems. There is a bulk of information available for certain values of q, namely 
when q — 1,2 and q is sufficiently large. In addition, the case d = 2 is special, and 
we shall return to this in Section 6. We assume throughout this section that q >1. 

Little is known in general about the numerical values of Pc{q)- For example, 
it is known that Pc{q) is Lipschitz-continuous and strictly increasing when d > 2, 
[68], and there is a striking conjecture (OP) that Pc{q) = \fql (1 + \fq) when d = 2 
(see Section 6.2). Some concrete inequalities involving the Pc(?) are implied by the 
comparison inequalities of Theorem 3.2. 

5.1 The subcritical phase, p < pdq) 

As remarked in Section 4.2, there is a unique random-cluster measure when p < 
Pc{q), and we shall denote this by 4>p,q- 

The key theorem for understanding the subcritical phase of percolation states 
that long-range connections have exponentially decaying probabilities. Such a result 
is believed to hold for all random-cluster models with q > 1, but no proof has been 
found (OP) which is valid for all > 1 and all p < Pc{q)- The full result is known 
only when g = 1, g = 2, or g is sufficiently large, and the three sets of arguments 
for these cases are somewhat different from one another. As for results valid for 
all q (> 1), the best that is currently known is that the connectivity function 
decays exponentially for sufficiently small p (this follows by Theorem 3.2 and the 
corresponding q = 1 result), and that it decays exponentially whenever it decays at 
a sufficient polynomial rate. We describe the last result next. 

As a preliminary we introduce another definition of a critical point. Let B{n) 
be the cube [— n,n]'^. We write 

Y{p,q) = limsup{n^-Vp,g(0 ^ dB{n))\ 

and Pg{q) = sup{p : Y{p, q) < ooj. Evidently < Pg{q) < Pc{q)i and it is believed 
that pg(g) = pc(g) for all g > 1 (OP). 
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Theorem 5.1 [76]. Let q>l,d>2, and < p < Pg{q)- There exists 7 = 7(p, q) 
satisfying 7 > such that 

4'p,q{0 ^ dB{n)) < e~'^'^ for all large n. 

The spirit of the theorem is close to that of Hammersley [84] and Simon-Lieb 
[112, 136] who proved exponential estimates when q = 1,2 subject to a hypothesis 
of finite susceptibility (that is, under the hypothesis that ^a;</'p,g(0 ^ x) < 00). 
The latter assumption is slightly stronger than the assumption of the above theorem 
when d = 2. 

Connectivity functions are expected to decay exponentially with a correction 
term of power order. More specifically, it is expected as reported in Section 4.2 
that 

(/>p,g(0 ^ a;) ~ ^^^3jy^exp(-|x|/0 as |a;| 00, 

for constants c{p,q) and ^{p,q), and for some suitable norm | • | on Z^. Such 
'Ornstein-Zernike' decay is a characteristic of many systems in their disordered 
phases. No proof is known (OP), except in the special cases when q = 1 and q — 2, 
[32, 33]. In [9] may be found a weaker result which bounds the fluctuations by a 
power-law when d = 2, under the assumption that the function does indeed decay 
exponentially. 

5.2 The supercritical phase, p > Pc{q) 

We assume as usual that g > 1, and we begin with a discussion of the number 
of infinite clusters. For uj e f2, let I{uj) be the number of infinite open clusters. 
Suppose that (^p^q is a translation- invariant member of 7?-p,g. If in addition (j)p^q is 
ergodic, then, by a well known theorem of Burton and Keane [30] , 

either 4>p,q{I = 0) = 1 or 4>p,q{I = 1) = 1; 

that is to say, the infinite open cluster is almost surely unique whenever it exists. It 
is noted in [30] that methods of ergodic decomposition enable the extension of such 
results to translation-invariant measures which are not necessarily ergodic. That 
is, under the assumption of translation-invariance alone, 

(/.j,,q(/e{o,i}) = i, 

which is to say that translation-invariant random-cluster measures have the 0/1- 
infinite-cluster property. A further comment on the use of ergodic decomposition 
in this context is to be found in [31]. 

In two dimensions, the supercritical process is best studied via the subcritical 
process which arises as its graphical dual (see Section 6). There are two general 
approaches to the supercritical phase in a general number d (> 3) of dimensions. 
The less powerful is to derive results for large p by comparison with percolation, 
the theory of which is relatively complete. Without an extra ingredient, such an 
approach will not reveal the structure of the supercritical phase all the way down 
to the critical value Pc{q)- As an example, we present one theorem concerning the 
uniqueness of random-cluster measures. 
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Theorem 5.2 [69]. If d > 2 and q > I, there exists p' = p'{d,q) < 1 such that 
^p,q ~ ^p,Q "whenever p > p' . 

It is an important open problem to prove that 0°^^ = 0^ ^ for aU p > Pc{q), 
or equivalently that there exists a unique random-cluster measure throughout the 
phase (OP). 

A more powerful approach, sometimes used in conjunction with the comparison 
argument summarised above, is the 'block argument' laid out in [36, 125]. One may 
think of block arguments as a form of rigorous renormalisation. One divides space 
into blocks, constructs events of an appropriate nature on such blocks, having large 
probabilities, and then allows these events to combine across space. There have 
been substantial successes using this technique, of which the most striking is the 
resolution, subject to certain side conditions, of the so-called Wulff construction for 
the asymptotic shape of large Ising droplets. 

Rather than discussing the physical background of the Wulff construction, we 
mention instead its impact on random-cluster models. Let B{n) = [—n,n]'^, and 
consider the wired random-cluster measure ^^^^-j ^ ^ with p > Pc{q)- The larger is 
an open cluster, the more likely it is to be joined to the boundary dB{n). Suppose 
that we condition on the event that there exists in B{'n) an open cluster C which 
does not touch dB{n) and which has volume of the order of the volume n'^ of the 
box. What can be said about the shape of C? Since p > Pc{q), there is little cost in 
having large volume, and the price is spent around its boundary. Indeed, the price 
may be expressed as a surface integral of an appropriate function termed 'surface 
tension'. This 'surface tension' may be specified as the exponential rate of decay 
of a certain probability. The Wulff prediction for the shape of C is that, when 
re-scaled in the limit of large n, its shape converges to the solution of a certain 
variational problem, that is, the limit shape is obtained by minimising a certain 
surface integral subject to a condition on its volume. 

No proof of this general picture for random-cluster models has appeared in the 
literature, although it is believed that the methods of [36, 37, 125] enable such 
a proof. The authors of [36] have instead concentrated on using random-cluster 
technology to solve the corresponding question for the asymptotic shape of large 
droplets in the Ising model. The outcome is an important 'large deviation' theorem 
which utilises block arguments and yields a full solution to the Ising problem when- 
ever the corresponding random-cluster model (which has q = 2) has parameter p 
satisfying p > Pc{2) and 0^ 2 — ^p,2- Here, Pc{2) is the limit of a certain decreas- 
ing sequence of critical points defined on slabs in Z^, and is conjectured (OP) to 
be equal to the critical point Pc{2). [Closely related results have been obtained 
in [24]. Fluctuations in droplet shape for two-dimensional random-cluster models 
have been studied in [10, 11].] 

The 'slab critical point' Pc{q) may be defined for any random-cluster model as 
follows. Fix Q > 1, and let d > 3. Let S{n,L) = [-n,n]'^-^ x [-L,L]. Let i;^;^ 
be the random-cluster measure on S{n, L) with parameters p, q (and with free 
boundary conditions). We denote by n(p, L) the property that: 

there exists a > such that, for all x e S{n, L) and all n, V'p,'^(0 ^ x) > a. 
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It is not hard to see that n(p, L) =^ n(p', L') Hp <p' and L < L' . It is thus natural 
to define the quantities 

Pc(g, L) = inf {p : n(p, L) occurs}, Pc{q) = hm pc(g, L), 
and it is clear that pdo) > PciQ)- 

Conjecture 5.3 [125]. Let q>l and d>3. We have that Pc{q) —Pcil)- 

Subject to a verification of this conjecture, and of a positive answer to the ques- 
tion of the uniqueness of random-cluster measures when p > pdq), the block ar- 
guments of [36, 125] may be expected to result in a fairly complete picture of the 
supercritical phase of random-cluster models with q > 1; see [37] also. 

The case = 1 is special, percolation enjoys a spatial independence not shared 
with general random-cluster models. This additional property has been used in 
the formulation of a type of 'dynamic renormalisation', which has in turn yielded a 
proof that Pc(l) — Pc{^) for percolation in three or more dimensions, [71, Chapter 
7, 74]. Such arguments do not to date have a random-cluster counterpart. 

As a further application of a block argument we note the following bound, [125], 
for the tail of the size of the open cluster C at the origin, 

</,J^^(|C| = n) < exp(-an(^-^)/^) for aU n, 

for some a = a{p,q) > 0, and valid for d > 3, 6 = 0, 1, and p sufficiently close to 1. 
The complementary inequality 

(Pl^g{\C\ =n)> exp(-a'n('^-^)/'^) for aU n, 

may be obtained for large p as done in the case of percolation, [71, Section 8.6]. 
5.3 Near the critical point, p ~ Pc{q) 

Surprisingly little is known about random-cluster measures near the critical point, 
except in the cases q = 1,2 and q large. In each such case, there are special 
arguments which are apparently not suitable for generalisation. We summarise 
such results as follows. 

I. Percolation, q = 1- 

There is a full theory of the subcritical and supercritical phases of percolation, 
[71]. The behaviour when p ~ ^^(1) has been the subject of deep study, and many 
beautiful results are known. Nevertheless, the picture is incomplete. For example, 
it is believed but not proved that 9{pc{l),l) = for all d > 2, but this is known 
only when d = 2 (because of special properties of two dimensions explored for 
in Section 6) and when d is large {d > 19 suffices) using a method termed the 'lace 
expansion'. The lace expansion explains also the values of some critical exponents 
when d is large; see, for example, [85, 86]. 

Great progress has been made in recent years towards understanding the phase 
transition when d = 2. The idea is to work at the critical point p = pc{l), and to 
observe the process over an increasing sequence of regions of . It is believed that 
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the process, re-scaled as the regions become larger, converges in a certain manner 
to a stochastic process generated in a prescribed way by a differential equation, 
known as a Lowner equation, which is driven in a certain way by a Brownian 
motion. Stochastic processes which arise in this way have been termed stochastic 
Lowner evolutions by Schramm, [135], and denoted SLE«, where k is the variance 
parameter of the Brownian motion. It is believed that the space of stochastic 
Lowner evolutions is a canonical family of processes which arise as scaling limits of 
discrete processes such as critical percolation, critical random-cluster models with 
Q < 4, self-avoiding walks, loop-erased random walk, and uniform spanning trees. 
Full proofs are not yet known (OP). We expand on this very important development 
in Section 6.4 

II. Ising model, q = 2. 

Integer values of q are special, and the value q = 2 particularly so because of certain 
transformations which permit the passage to a model which might be termed a 
'Poisson graph'. Let G = {V, E) be a finite graph and let < A < oo. Suppose that 
TT = {7r(e) : e G E} is a family of independent random variables each having the 
Poisson distribution with parameter A. We now construct a random graph = 
(V, Et^) having vertex set V and, for each e ^ E, having exactly 7r(e) edges in parallel 
joining the endvertices of the edge e [the original edge e is itself removed] . We call 
a Poisson graph with intensity A, and write and E;^ for the appropriate 
probability measure and expectation operator. 

We introduce next the concept of a flow on an oriented graph. Let g G {2,3,...} 
and let G' = {V , E') be a finite oriented graph. Let / : ^ {0, 1, 2, . . . , g - 1}. 
For X e V, the total flow into x is the sum of ±f{e') over all edges e' incident to x, 
with when e' is oriented towards x and —1 otherwise. The function / is called 
a mod-q flow if the total fiow into x is zero (modulo q) for all x & V. The mod-g 
flow / is called non-zero if f{e') ^ for every e' e E' . We write Fq{G') for the 
number of non-zero mod-g flows on G'. It is a remarkable fact, [142], that Fq{G') 
does not depend on the orientations of edges in E' ^ and thus one may define Fq(G') 
unambiguously for any unoriented graph G' . 

We return now to the Poisson graph Gj^. For x,y E V, x y, we denote by 
G^'y the graph obtained from G-,^ by adding an edge with endvertices x, y. [If x 
and y are already adjacent in Gtt, we add exactly one further edge between them.] 
Connection probabilities and flows are related by the following theorem, which may 
be proved using properties of Tutte polynomials (see [142] and Section 2.5). 

Let G = {V, E) be a finite graph, and write 4^G,p,q for the random-cluster measure 
on G with parameters p, q. 



Theorem 5.4 [63, 73]. Let q e {2,3, ...} and < p = 1 - e-^« < 1. We have 



This formula takes on an especially simple form when q — 2, since non-zero 
mod-2 flows necessarily take only the value 1. It follows that, for any graph G', 
^\{F2{G')) equals the PA-probability that the degree of every vertex of G' is even, 
[1]. Observations of this sort have led when q = 2 to the so-called 'random-current' 



that 



{q - l)(t)G,p,q{x ^ 



y) 



¥.^{Fq{G%^y)) 
^x{Fq{G^)) 



for all x,y & V, x ^ y. 
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expansion for Ising models, thereby after some work [1, 2, 5] leading to proofs 
amongst other things of the following, expressed here in the language of random- 
cluster measures. 

(i) When g = 2 and p < Pc{q), we have exponential decay of the radius distribu- 
tion, 

(/>p,2(rad(C) = n) < e"''^ for aU n, 

where rj — r](j)) > 0; exponential decay of the two-point connectivity function 
follows. 

(ii) When q = 2 and d ^ 3, there is a unique random-cluster measure 0^,2 for all 
p, in that |7^p,q| = 1. 

(iii) The phase transition is continuous when q = 2 and d ^ 3. In particular, 
6'^(pc(2), 2) = 6'^(pc(2),2) = 0, and the edge-density h{p,2) = 0p,2(e is open) 
is a continuous function of p at the critical point Pc{2). 

(iv) When d > 4, some (at least) critical exponents take their mean-field values, 
and depend no further on the value of d. 

Note that the nature of the phase transition in three dimensions remains curi- 
ously undecided (OP). 

III. The case of large q. 

It is not known whether the phase transition is continuous for all small q (OP). 
The situation for large q is much better understood owing to a method known as 
Pirogov-Sinai theory [123, 124] which may be adapted in a convenient manner to 
random-cluster measures. The required computation, which may be found in [105], 
has its roots in an earlier paper [103] dealing with Potts models. A feature of such 
arguments is that they are valid 'all the way to the critical point' (rather than for 
'small p' or 'large only), so long as q is sufficiently large. One obtains thereby a 
variety of conclusions including the following. 

(i) The edge-densities h^{p, q) — (pp^q{e is open), 6 = 0, 1, are discontinuous func- 
tions of p at the critical point. 

(ii) The percolation probabilities satisfy 0^{pc{q): q) — 0, 9^{pc{q), q) > 0. 

(iii) There is a multiplicity of random-cluster measures when p = Pc{q), in that 

(iv) If p < Pc{q), there is exponential decay and a mass gap, in that the unique 
random-cluster measure satisfies 

0p,g(O ^ en) = e-(^+°(^))"/^ as n ^ oo, 

where e„ = (n, 0, 0, . . . , 0) and the correlation length ^ = ^{p, q) is such that 

1™pTpc(<j) ^(P: Q) = ^(q) < oo- 

(v) If d — 3 and p — Pciq), there exists a non-translation-invariant random-cluster 
measure, [38, 116]. 

It is not especially fruitful to seek numerical estimates on the required size Q{d) 
of q for the above conclusions to be valid. Such estimates may be computed, but 
turn out to be fairly distant from those anticipated, namely Q{2) = A, Q{d) = 2 for 
d>6. 
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The proofs of the above facts are rather comphcated and will not be explained 
here. Proofs are much easier and not entirely dissimilar when c/ = 2, and a very 
short sketch of such a proof is provided in Section 6.3. 



6. In two dimensions 

The duality theory of planar graphs provides a technique for studying random- 
cluster models in two dimensions. We shall see in Section 6.1 that, for a dual pair 
(GjC^) of planar graphs, the measures 4>G,p,q ci-nd 4>G'^^p<^^q are dual measures in 
a certain geometrical sense, where p, p"^ are related by — p'^) = q(1 — p)/p- 

Such a duality permits an analysis by which many results for may be derived. 
Of particular interest is the value of p for which p = p^. This 'self-dual point' is 
easily found to be p = Psd(Q') where 

Psd(g) = 



and it is conjectured that Pc{(l) = Padiji) for g > 1. 
6.1 Graphical duality 

Let G = (V, E) be a simple planar graph imbedded in M^. We obtain its dual graph 
= {V^j E"^) as follows (the roman letter 'd' denotes 'dual' rather than number 
of dimensions). We place a dual vertex within each face of G, including the infinite 
face of G if G is finite. For each e E E we place a dual edge — {x'^,y'^) joining 
the two dual vertices lying in the two faces of G abutting e; if these two faces are 
the same, then x'^ = y'^ and is a loop. Thus E'^ is in one-one correspondence to 
E. It is easy to see that the dual of is isomorphic to L^. What is the relevance 
of graphical duality to random-cluster measures on G? 

Suppose that G is finite. Any configuration G O (= {0, 1}^) gives rise to a dual 
configuration lo^ lying in the space — {0, l}^"* defined by a;'^(e'^) = 1 — uj{e). As 
before, to each configuration oj'^ corresponds the set r]{uj^) — {e*^ e E'^ : LO^{e'^) = 
1} of its 'open edges'. Let f{uj) be the number of faces of the graph (F, 77(0;)), 
including the infinite face By drawing a picture, one may easily be convinced (see 
Figure 6.1) that the faces of {V,r]{uj)) are in one-one correspondence with the 
components of {V^ , r]{u!'^)) , and therefore f{u) — k{u>'^), in the obvious notation. 
We shall make use of Euler's formula (see [147]), 

k{u;) = \V\-\r]{uj)\+ f{u;) -1, ou e n. 
The random-cluster measure on G is given by 

Using Euler's formula and the equality /(cu) = k{uj^), we find that 

0G,p,g(^) = 0Gd,pd,g(a;'^) for oj e fi, 
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Figure 6.1. A primal configuration uj (with solid lines and vertices) and its dual 
configuration u)'^ (with dashed lines and hollow vertices). The arrows join the given 
vertices of the dual to a dual vertex in the infinite face. Note that each face of the 
primal graph (including the 'infinite face') corresponds to a unique component of the 
dual graph. 

where the dual parameter is given according to 

1 — p"^ p 

The unique fixed point of the mapping p i— > p*^ is given by p = Psd(?) where 
Psdil) — ■\/q/ (1 + v^)- note at this point that 

by Euler's formula. This representation of the random-cluster measure at the 'self- 
dual point' Psd{(l) highlights the duality of measures. 

Turning to the square lattice, let = [0,n]^, whose dual graph AJ^ may be 
obtained from [— 1, n]^ + (|, ^) by identifying all boundary vertices. By the above, 

for configurations lo on A^ (and with a small 'fix' on the boundary of A^). Letting 
n — > oo, we obtain that ^(A) — (p^^ qi^'^) ^r all cylinder events A, where 
A"^ = {u'^ : u; e A}. 

6.2 Value of the critical point 

Consider the random-cluster process on the two-dimensional lattice = (Z^,E^), 
with parameters p and q satisfying q > 1. The following remarkable conjecture is 
widely believed (OP). 

Conjecture 6.1. Let q > 1. The critical value Pc{q) o/L^ is given by 
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This conjecture is known to hold when q = 1, q = 2, and for q > 25.72. The q = 1 
case was answered by Kesten [97] in his famous proof that the critical probability of 
bond percolation on is ^. For q = 2, the value of Pc{2) given above agrees with 
the celebrated calculation by Onsager [119] of the critical temperature of the Ising 
model on , and is implied by probabilistic results in the modern vernacular of [2] . 
The formula for Pc{q) has been established rigorously in [104, 105] for sufficiently 
large (real) values of q, specifically q > 25.72 (see also [70]). 

Conjecture 6.1 arises in a natural manner from the observation that is a self- 
dual graph, and p — Psdil) — \/?/ (1 + ^/q) is the self-dual point of a random-cluster 
measure on with parameters p, q. 

Several other remarkable conjectures about the phase transition in may be 
found in the physics literature (see [14]), as consequences of 'exact' but non-rigorous 
arguments involving ice-type models. These include exact formulae for the as- 
ymptotic behaviour of the partition function limA|z2{-^A.p,g}^^''^', and also for the 
edge-densities h^{Psd{Q),Q) — 0p ^(g) is open), b — 0,1, at the self-dual point 
Psd(g)- 

Progress towards a rigorous verification of the conjecture may be summarised 
briefly as follows. Using an argument, [152], taken from percolation using the 
uniqueness of infinite open clusters, we obtain by duality that 6^{psd{Q),Q) = 
(see [69, 143]), whence the critical value of the square lattice satisfies Pc{q) > Psd{<l) 
for q > 1. The complementary inequality Pc(?) ^ Psd{Q) has eluded mathematicians 
despite progress by physicists, [87]. 

Suppose on the contrary that Pc{q) > Psdil), so that pdl)'^ < Psdil)- For 
P G iPcil)'^ T Pcig)) we have also that p^ e iPciQ)^ ■, Pc{<l)) ■ Therefore, for p e 
iPcil)'^ T Pci^)) , both primal and dual processes comprise, almost surely, the union of 
finite open clusters. This contradicts the intuitive picture, supported for p ^ pdq) 
by our knowledge of percolation, of finite clusters of one process fioating in an 
infinite ocean of the other process. 

Exact values for the critical points of the triangular and hexagonal lattices may 
be conjectured similarly, using graphical duality together with the star-triangle 
transformation, [14, 101]. 

Rigorous numerical upper bounds of impressive accuracy have been achieved for 
the square lattice and other two-dimensional lattices via an extension of the basic 
model to a larger class termed in [8] the 'asymmetric random-cluster model'. The 
bound in question for is 

^/l-q ^ + ^ 

For example, when q = 10, we have that 0.760 < Pc(lO) < 0.769, to be compared 
with the conjecture that Pc(lO) = \/T0/ (1 + vTO) ~ 0.760. A valuable consequence 
of the comparison methods developed in [8] is the exponential decay of connectivity 
functions when q > 2 and p is such that 

p < Psdiq - 1) = 
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6.3 First-order phase transition 

There is a special argument discovered first for Potts models, [104], which may be 
used to show first-order phase transition when q is sufficiently large. 

Let a,j be the number of self-avoiding walks on beginning at the origin. It 
is standard, [114], that a^"^ — as n — > cxd, for some constant n called the 

connective constant of the lattice. Let Q = { ^ (a* + a//U^ — 4) } . We have that 
2.620 < /X < 2.696 (see [137]), whence 21.61 <Q < 25.72. We set 

1, / (i+y/gr i 

^^^) = 24^°n^;^|, 

noting that '^{q) > if and only \i q> Q. We write B{ri) = [—n, n]^. 

Theorem 6.2 [70, 104]. Ifd — 2 and q > Q then the following hold. 

(a) The critical point is given by Pc{q) = ■\/q/ (1 + -s/q)- 

(b) We have that 9^{pc{q),q) > 0. 

(c) For any ^/j < il^{q) and all large n, (t^'^^i^q-^ ^ dB{n)) < e""''^. Hence, in 
particular, 0'^{Pc{q),q) — 0. 

The idea of the proof is as follows. There is a partial order on circuits of 
given by r < r' if the bounded component of \ F is a subset of the bounded 
component of \ F'. We work at the self-dual point p — Psdil), and with the box 
B{n) with wired boundary conditions. An 'outer contour' is defined to be a circuit 
F of the dual graph B{n)'^ all of whose edges are open in the dual (that is, they 
traverse closed edges in the primal graph B{n)), and which is maximal with this 
property. Using self-duality, one may show that 

0B(n),p3d(5),g(r is an outer circuit) < ^ (^J^^^T^ 

for any given circuit F of B{n)'^. Combined with a circuit-counting argument of 
Peierls-type involving the connective constant, this estimate implies after a little 
work the claims of Theorem 6.2. The idea of the proof appeared in [104] in the 
context of Potts models, and the random-cluster formulation may be found in [70] . 

We stress that corresponding conclusions may be obtained for general d {> 2) 
when q is sufficiently large {q > Q{d) for suitable Q{d)), as shown in [105] using 
so-called Pirogov-Sinai theory. Whereas, in the case d = 2, the above duality 
provides an especially simple proof, the proof for general d utilises nested sequences 
of surfaces of M'^ and requires a control of the effective boundary conditions within 
the surfaces. 

6.4 SLE LIMIT WHEN g < 4 

Many exact calculations are 'known' for critical processes in two dimensions, but the 
physical arguments involved have sometimes appeared in varying degrees magical or 
revelationary to mathematicians. The new technology of stochastic Lowner evolu- 
tions (SLE), discovered by Schramm [135] and mentioned in Section 5.3, threatens 
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to provide a rigorous underpinning of many such arguments in a manner most 
consonant with modern probabihty theory. Roughly speaking, the theory of SLE 
informs us of the correct weak Hmit of a critical process in the limit of large spatial 
scales, and in addition provides a mechanism for performing calculations for the 
limit process. 

Let H = (—00,00) X (0,00) be the upper half-plane of IR^, with closure H. We 
view H and H as subsets of the complex plane. Consider the ordinary differential 
equation 

dt gt{z) - B^t 

subject to the boundary condition go{z) = z, where t e [0, 00), k is a positive 
constant, and {Bf : t > 0) is a standard Brownian motion. The solution exists 
when gt{z) is bounded away from B^f More specifically, for 2; G H, let be the 
infimum of all times r such that is a limit point of gsiz) — B^g in the limit as 
s t T. We let 

Ht = {z en:T^>t}, Kt = {zeW:T^< t}, 

so that Hi is open, and is compact. It may now be seen that gt is a conformal 
homeomorphism from Ht to H. 

We call ((7t : t > 0) a stochastic Lowner evolution (SLE) with parameter k, 
written SLE^, and we call the Kt the hulls of the process. There is good reason 
to believe that the family K — {Kt : t > 0) provides the correct scaling limit of a 
variety of random spatial processes, the value of k, being chosen according to the 
process in question. General properties of SLE^j, viewed as a function of k have 
been studied in [130, 145], and a beautiful theory has emerged. For example, the 
hulls K form almost surely a simple path if and only ii k < 4. If ac > 8, then SLE^j 
generates almost surely a space-filling curve. 

Schramm [135] has identified the relevant value of n for several different processes, 
and has indicated that percolation has scaling limit SLEg, but full rigorous proofs 
are incomplete. In the case of percolation, Smirnov [138, 139] has proved the very 
remarkable result that, for site percolation on the triangular lattice, the scaling 
limit exists and is SLEe (this last statement is illustrated and partly explained in 
Figure 6.2), but the existence of the limit is open for other lattices and for bond 
percolation. 

It is possible to perform calculations on stochastic Lowner evolutions, and in 
particular to confirm, [110, 140], the values of many critical exponents associated 
with percolation (or, at least, site percolation on the triangular lattice). The con- 
sequences are in agreement with predictions of mathematical physicists previously 
considered near- miraculous (see [71], Chapter 9). In addition, SLEg satisfies the 
appropriate version of Cardy's formula, [34, 107]. 

The technology of SLE is a major piece of contemporary mathematics which 
promises to explain phase transitions in an important class of two-dimensional dis- 
ordered systems, and to help bridge the gap between probability theory and confor- 
mal field theory. It has already provided complete explanations of conjectures, by 
mathematicians and physicists, associated with two-dimensional Brownian motions 
and specifically their intersection exponents and fractionality of frontier, [108, 109]. 
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Figure 6.2. Site percolation on the triangular lattice with p equal to the critical 
point ^, and with a mixed boundary condition along the lower side. The interface 
traces the boundary between the white and the black clusters touching the boundary, 
and is termed the 'exploration process'. In the limit of small lattice-spacing, the 
interface converges in a certain manner to the graph of a function which satisfies the 
Lowner differential equation driven by a Brownian motion with variance parameter 
K= Q. 

Extra work is needed in order to prove the validity of the hmiting operation 
for other percolation models and random processes. In another remarkable recent 
paper [111], Lawler, Schramm, and Werner have verified the existence of the scaling 
limit for loop-erased random walk and for the uniform spanning tree Peano curve, 
and have shown them to be SLE2 and SLEg respectively. It is believed that self- 
avoiding walk on L^, [114], has scaling limit SLEg/s. 

We turn now to the random-cluster model on L'^ with parameters p and q. 
For 1 < < 4, it is believed that the percolation probability 0{p,q), viewed as 
a function of p, is continuous at the critical point Pciq) (OP)? and furthermore 
that Pc{q) = \^/{^ + v^)- seems likely that, when re-scaled in the manner 
similar to that of percolation (illustrated in Figure 6.2), the exploration process 
of the model converges to a limit process of SLE type. It then remains only to 
specify the parameter k of the limit in terms of q. It has been conjectured in 
[130] that K satisfies cos(47r/K) = — |y^, k G (4,8). This value is consistent with 
Smirnov's theorem [139], and also with the finding of [111] that the scaling limit 
of the uniform spanning tree Peano curve is SLEg, on recalling that the uniform 
spanning tree measure is obtainable as a limit of the random-cluster measure as 
p,q 10 (recall Section 2.4). 

There are uncertainties over how this programme will develop. For a start, the 
theory of random-cluster models is not as complete as that of percolation and of 
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the uniform spanning tree. Secondly, the existence of spatial limits is currently 
known only in certain special cases. The programme is however ambitious and full 
of promise, and should ultimately yield a full picture of the critical behaviour — 
including values of exponents — of random-cluster models, and hence of Ising/Potts 
models, with q < 4:. 



While considerations of 'real space-time' support the study of such models on 
lattices such as L'^, it has proved rewarding also to analyse the random-cluster model 
on certain other graphs including complete graphs and trees. It is the presence of 
circuits in the underlying graph which is the root cause of dependence between 
the states of edges, and for this reason it is the complete graph which provides 
an appropriate setting for what is termed 'mean-field theory', in which vertices 
'interact' with all other vertices rather than with a selected subset of 'neighbours'. 
Trees, on the other hand, contain no circuits, and their random-cluster theory is 
therefore sterile unless one introduces boundary conditions. [A different approach 
to mean- field theory has been studied in [99], namely on L'^ for large d.] 

7.1 On complete graphs 

The mean-field Potts model may be formulated as a Potts model on the complete 
graph Kji, being the graph with n labelled vertices every pair of which is joined by 
an edge. The study of such a process dates back at least to 1954, [100], and has 
been continued over the last fifty years [26, 99, 151]. The model is exactly soluble 
in the sense that quantities of interest may be calculated exactly and rigorously. It 
is therefore not surprising that the corresponding random-cluster models (for real 
q) have 'exact solutions' also, [26]. 

Consider the random-cluster measure il^n,\,q = 4'Kn,x/n,q on the complete graph 
Kn, having parameters p — X/n and q; this is the appropriate scaling to allow an 
interesting limit as n ^ oo. In the case q — I, this measure is product measure, 
and therefore the ensuing graph is an Erdos-Renyi random graph [25, 90]. The 
overall picture for general values of q is rather richer than for the case q = 1, and 
many exact calculations may be performed rigorously. It turns out that the phase 
transition is of first-order if and only if q > 2, and the behaviour of the system 
depends on how A compares with a 'critical value' Xc{q) taking the value 



From the detailed picture described in [26] the following information may be 
extracted. The given properties occur with 7/;„^A,g-probability tending to 1 as n — > 

CO. 

I. Subcritical case, when X < Xc{q)- 

The largest component of the graph is of order logn. 

II. Supercritical case, when X > Xc{q). 



7. On complete graphs and trees 
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Figure 7.1. The function 6{X, q) for the three cases q < 2, q = 



2,q>2. 



There is a 'giant component' having order ^(A, q)n where 9 is defined to be the 
largest root of the equation 



IV. Critical case, when X — Xc{q), q > 2. 

The largest component is either of order logn or of order ^(A, q)n, where 9 is given 
as in case II above. 

The dichotomy between first- and second-order phase transition is seen by study- 
ing the function 0(A, q), sketched in Figure 7.1. When < q < 2, the function 0(A, q) 
descends continuously to as A J, Xc{q). On the other hand, this limit is strictly 
positive when q > 2. 

The above results are obtained via a relationship between the model for general 
q and the model for the special case q = I. The latter system has been analysed 
extensively, [25, 90]. We illustrate the argument in the case q > 1; a similar 
approach is valid when q < 1. Consider the open clusters Ci, C2, ■ ■ ■ , Cm of a sample 
from the random-cluster measure 0K„,p,q- We colour each such cluster red with 
probability p, and white otherwise, different clusters receiving independent colours. 
We delete all vertices in white clusters, and let H denote the remaining graph, 
comprising a certain random number N of vertices (from the red clusters) together 
with certain open edges joining pairs of them. It may be seen that, conditional 
on the value of A^, the measure governing H is the random-cluster measure with 
parameters p and qp. We choose p = 1/q to obtain an Erdos-Renyi random graph 
on a random set of vertices. This is the observation which permits the full analysis 
to proceed. 

One consequence of this study is an explicit identification of the exponential 
asymptotics of the partition function Zx„,\/n,qj namely of the limit function 



This provides information via the Gartner-Ellis theorem, [43] , concerning the large- 
deviation theory of the number of clusters in such systems. 



l + {q-l)9 
1-9 



III. Critical case, when X = Xdq) , < q < 2. 
The largest component has order n^/*^. 
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7.2 On trees and non- amenable graphs 

Whereas physical considerations support the study of interacting systems on finite- 
dimensional lattices, mathematicians have been drawn also to the study of general 
graphs, thus enabling a clearer elucidation of the mathematical structure of such 
systems (see the discussion in [134]). A subject of special focus has been the class 
of graphs for which the ratio of surface to volume of finite boxes does not approach 
zero in the limit as the size of the box tends to infinity. A prime example of such 
a graph is an infinite regular tree with vertex degree at least three. We make the 
distinction more concrete as follows. Let G = {V, E) be an infinite connected graph 
with finite vertex degrees. For C we define its boundary dW to be the set 
of all w e having some neighbour v not belonging to W . The countably infinite 
graph G is called amenable if its 'Cheeger constant' 



satisfies k.{G) = 0; G is called non- amenable if k.{G) > 0. It is easily seen that 
is amenable, whereas an infinite regular tree with degree at least three is non- 
amenable. 

The role of amenability in probability theory has been evident since the work 
of Kesten [95, 96] concerning random walks on a general graph G. More relevant 
to this review are [17, 75], which consider the number of infinite clusters in the 
bond percolation model on G. Suppose G is a quasi-transitive graph (that is, its 
vertex set has only finitely many orbits under its automorphism group). Suppose 
in addition that G is amenable. Consider bond percolation on G with density p. It 
may be proved as in [30, 58] that the number / of infinite open clusters satisfies 



That is, if an infinite open cluster exists, then it is almost surely unique. Under 
similar assumptions on a non-amenable graph G, it is believed but not yet proved 
in full generality that there exists an interval of values of p for which Pp{I = 
oo) = 1; see, for example, the discussion in [113]. A corresponding question for 
random-cluster models is to ascertain for which graphs G and values of q there is 
non-uniqueness of random-cluster measures for an interval of values of p. [Recall 
Theorem 3.7, easily extended to more general amenable graphs, which states that, 
for g > 1, there is a unique random-cluster measure on L'' for all except at most 
countably many values of p.] See [82, 92] and especially [134, Section 6.1] for recent 
accounts of this and associated questions, and [80] for an analysis of random-cluster 
measures on regular trees. 



Let be a probability measure on a space (n,jF). We may study stochastic pro- 
cesses taking values in Q which converge weakly to in the limit of large times. 
There are a multiplicity of reasons for and benefits in studying time-evolutions. 




either Pp{I = 0) = 1 or Pp{I = 1) = 1. 



8. Time-evolutions of random-cluster models 
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First, physical systems generally have dynamics as well as equilibria. Secondly, 

new questions of interest arise, such as that of the estimation of a relaxation time. 
Thirdly, the dynamics thus introduced can yield a new technique for studying the 
limit measure fx. 

When studying a physical system, it is often acceptable to restrict oneself to dy- 
namics which are reversible in time. In Section 8.1, we describe a natural reversible 
dynamic for a random-cluster model, akin to the Glauber dynamics of the Ising 
model. This dynamic permits an extension which couples together the random- 
cluster measures on a given graph as p and q range over their possible values. 

The problem commonly arises in statistics, computer science, and statistical 
physics of how to obtain a sample from a system governed by a probability measure 
with complex structure. In Section 8.2 we summarise the Propp-Wilson 'coupling 
from the past' approach, [128], to this problem in the context of the random-cluster 
measure. 

Since Potts models may be obtained from random-cluster models, there is an 
interplay between the dynamics for these two systems. A famous instance of this 
relationship is the so-called Swendsen-Wang dynamic [141], which is described in 
Section 8.3. 

We assume throughout this section that G = {V, E) is a, finite connected graph, 
and that (j)p^q is the random-cluster measure on n = {0, 1}^ with < p < 1 and 
q>0. 

8.1 Reversible dynamics 

We shall consider transitions from a configuration a; to configurations which differ 
from u) on one edge only. Thus we introduce the following notation. For a; G O 
and e E E, let and Ue be the configurations obtained by 'switching e on' and 
'switching e off', respectively, that is 



Let {Xt : t > 0) be a Markov chain, [77], on the state space fl with generator 
Q = {<lui,u}' : a;, a;' e Cl} satisfying 



where D{e, ^) is the indicator function of the event that the endpoints of e are joined 
by no open path of ^. This specifies the rate at which single edges are acquired or 
lost. We set q^^^^ = if c*; and ^ differ on two or more edges, and we choose the 
diagonal elements qu,,uj in such a way that Q, when viewed as a matrix, has row 
sums zero, that is. 




{l-p)q 



,^(e,u;e) 



for a; e ri, e e E, 




It is elementary that the 'detailed balance equations' 



u!,uj' e Q, 
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hold, whence X is reversible with respect to (pp^q. It follows by the irreducibility of 
the chain that Xf =^ (f>p^q as t ^ cxd (where denotes weak convergence). There 
are of course many Markov chains with generators satisfying the above detailed 
balance equations, the important quantity is the ratio quj,uj' /Quj',u)- 

Two extensions of this dynamical structure which have proved useful are as 
follows. The evolution may be specified in terms of a 'graphical representation' 
constructed via a family of independent Poisson processes. This allows a natural 
coupling of the measures 0p ^ for different p and q. Such couplings are monotone in 
p when g > 1. One may similarly couple the unconditional measure (t>p,q{-) and the 
conditioned measure (/>p^g(- | A). Such couplings permit probabilistic interpretations 
of differences of the form (f)pi^q{B \ A) — (j)p^q{B) when q > l, p < p' , and A and B 
are increasing, and this can be useful in particular calculations (see [19, 68, 69]). 

We turn now to the thermodynamic limit, and the question of the structure of 
a Markovian random-cluster process on an infinite connected graph. In the case 
q > 1, the above couplings are monotone in the choice of the underlying graph G. 
Therefore there exist 'limit dynamics' as G passes through an increasing sequence 
of finite graphs. Boundary conditions may be introduced, and one may obtain 
thereby a certain Markov process ^ = {(t '■ t > 0) on the state space [0, 1]^, where 
E is the limiting (infinite) edge set. This process, which does not generally have 
the Feller property, generates a pair of 'level-set processes' taking values in {0, l}'^, 
defined for < p < 1 by 

Cr"(e) = l{Ct(e)>i-p}, Cr^(e) = l{a(e)>i-p}, eeE, 

where, as before, 1a denotes the indicator function of an event A. These two 
processes arc Markovian and are reversible with respect to the infinite-volume free 
and wired random-cluster measures, respectively. See [69]. 

Note that the generator of the Markov chain given above depends on the random 
variable D(e,uJe), and that this random variable is a 'non-local' function of the 
configuration co in the sense that there is no absolute bound on the distance from e 
of edges whose states may be relevant to its value. It is this feature of non-locality 
which leads to interesting complications linked in part to the 0/1-infinite-cluster 
property introduced before Theorem 3.9. Further discussion may be found in [69, 
121]. 

8.2 Coupling from the past 

In running Monte Carlo experiments, one requires the ability to sample from the 
probability measure The Markov chain Xt of Section 8.1 certainly converges 
weakly to (f)p^q as t — > oo, but this is not as good as having in the hand a sample 
with the exact distribution. Random-cluster measures are well suited to the Propp- 
Wilson approach to sampling termed 'coupling from the past', [128], and we sketch 
this here. Some illustrations may be foimd in [91]. 

First we provide ourselves with a discrete-time reversible Markov chain {Z^ : 
n > 0) on the state space Q having invariant measure 4>p,q. The so-called heat-bath 
algorithm provides a suitable example of such a chain, and proceeds as follows. At 
each stage, we pick a random edge e, chosen uniformly from E and independently 
of all earlier choices, and we make e open with the correct conditional probability, 
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given the configuration on the other edges. The corresponding transition matrix is 
given by n = {tTj^^i^/ : a;,a;' e Q} where 

TT, — 



\E\ (/.p,,(a;-) + (/.p,,(a;e)' 

_ 1 0p,q(<^e) 



A neat way to do this is as foUows. Suppose that — u>. Let e„ be a random edge 
of E, and let Un be uniformly distributed on the interval [0, 1], these variables being 
chosen independently of all earlier choices. We obtain Z^+i from u by retaining 
the states of all edges except possibly that of e„. We set 

Zn+i{en) = if and only if Un < '^P''?^^^) 



Thus the evolution of the chain is determined by the sequences Cn, Un, and the 
initial state Zq. One may make this construction explicit by writing Zn+i = 
'ip{Zn,en,Un) (= ip{(jiJ , en,Un)) for some function ip : Q, x E x [0,1] ^ fl. It is 
easily seen by the HoUey condition of Section 3.1 that, if > 1, and for every e and 
u, the function e, u) is non-decreasing in its first argument. It follows that the 
coupling is 'monotone' in the sense that, if a; < a;', then the chain starting at u) lies 
at all times beneath the chain starting at ui' (using the partial order on Q). 

We let W = {W{ui) : a; e il) be a vector of random variables such that 
W{uj) has the distribution of Z\ conditional on Zq = uj. Following the scheme 
described above, we may take W{uj) = ifj{u,e,U) where e and U are chosen at 
random. Let W-m, m > 1, be independent random vectors distributed as W, 
that is, W-m{-) — '*/'(■; e^, C/^) where the set {(e^,C/^) : m > 1} comprises in- 
dependent pairs of independent random variables, each Cj being uniform on E, 
and each Ui being uniform on [0,1]. We now construct a sequence Y-n, > 1, 
of random maps from to O by the following inductive procedure. First, for 

G O, we set y_i(aj) = W-i{uj). Having found YLi, Y'_2, . . . , l^-m, we define 
= Y_ni{W-m-i{i^)) ■ That is, is obtained from u> by passing 

in one step to W-m-i{<jij), and then applying to this new state. The exact 
dependence structure of this scheme is an important ingredient of what follows. 

We stop this process at the earliest time m at which 'coalescence has occurred', 
that is, at the moment M given by M = min{m : Y-ni{-) is the constant function}. 
It is a theorem, [128], that M is (/)p^g-a.s. finite and, for any lv, the random output 
Y-Miuj) is governed exactly by the probability measure 4>p,q. 

This procedure looks unwieldy, since Cl is typically rather large, but the reality 
is simpler when q > 1. By the monotonicity of the above coupling when q > 1, 
it suffices to follow the trajectories of the 'smallest' and 'largest' configurations, 
namely those beginning, respectively, with every edge closed and with every edge 
open. The processes starting at intermediate configurations remain sandwiched 
between the extremal processes, for all times t. Thus one may define M by M = 
min{m : y_^(0) = y_^(l)}, where and 1 denote the vectors of zeros and ones as 
before. 
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8.3 Swendsen-Wang dynamics 

It is a major target of statistical physics to understand the time-evolution of dis- 
ordered systems, and a prime example lies in the study of the Ising model. A 
multiplicity of types of dynamics have been proposed. The majority of these share 
a quality of 'locality' in the sense that the evolution involves changes to the states 
of vertices in close proximity to one another, perhaps single spin-flips, or spin- 
exchanges. The state space is generally large, of size 2^ where N is the number 
of vertices, and the Hamiltonian has complicated structure. When subjected to 
'local dynamics', the process may approach equilibrium very slowly (see [115, 133] 
for accounts of recent work of relevance). 'Non-local dynamics', on the other hand, 
have the potential to approach equilibrium faster, since they permit large jumps 
around the state space, relatively unconstrained by neighbourly relations. The 
random-cluster model has played a role in the development of a simple but attrac- 
tive such system, namely that proposed by Swendsen and Wang [141] and described 
as follows for the Potts model with q states. 

As usual, G — {V, E) is a finite graph, typically a large box in Z*^, and E = 
{1,2,..., q}^ is the state space of a Potts model on G. We write Q, = {0, 1}^. 
Suppose that, at some time n, we have obtained a configuration an (e E). We 
construct cr„+i as follows. Let p=l — e~^^ where < f3J < oo. 

I. We let a;„ e f2 be given as follows. For e = {x, y) e 



different edges receiving independent states. The edge configuration ojn is 
carried forward to the next stage. 

II. To each cluster C of the graph {V, rj{ujn)) we assign an integer chosen uniformly 
at random from the set {1, 2, . . . , q'}, different clusters receiving independent 
labels. We let a"„_|_i(a;) be the value thus assigned to the cluster containing 
the vertex x. 

It may be checked that the Markov chain {an : n > 0) has as unique invariant 
measure the Potts measure on E with parameters P and J. (Recall paragraph (c) 
of Section 2.3.) 

The Swendsen-Wang algorithm leads to samples which generally converge to 
equilibrium faster than those defined via local dynamics. This is especially evident 
in the 'high (3' (or 'low temperature') phase, for the following reason. Consider 
for example the simulation of an Ising model on a finite box with free boundary 
conditions, and suppose that the initial state is -|-1 at all vertices. If /? is large, 
then local dynamics result in samples which remain close to the '+ phase' for a very 
long time. Only after a long wait will the process achieve an average magnetisation 
close to 0. Swendsen-Wang dynamics, on the other hand, can achieve large jumps 
in average magnetisation even in a single step, since the spin allocated to a given 
large cluster of the corresponding random-cluster model is equally likely to be either 
of the two possibilities. A rigorous analysis of rates of convergence is however 



if an{x) ^ an{y), let u)n{e) = 0, 




(c) Geoffrey Grimmett, 1 February 2008 



40 



GEOFFREY GRIMMETT 



incomplete. It turns out that, at the critical point, Swendsen-Wang dynamics 
approach equilibrium only slowly, [28]. A further discussion is available in [61]. 

Algorithms of Swendsen-Wang type have been described for other statistical 
mechanical models having graphical representations of random-cluster-type; see 
[41, 42]. Related work may be found in [149]. 
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